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ABSTRACT 

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the 
missing baryons. However, the kSZ signal is overwhelmed by various contaminations and 
the cosmological application is hampered by loss of redshift information due to the projec- 
tion effect. We propose a kSZ tomography method to alleviate these problems, with the aid 
of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar ve- 
locity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the 
product projected along the line of sight with a proper weighting as an estimator of the true 
kSZ temperature fluctuation 0. Since the underlying directional dependence in the estimator 
closely resembles that in the true kSZ signal 0, is tightly correlated with 0. It thus avoids 
the problem of null correlation between the galaxy density and 0, which prohibits the kSZ 
extraction through the usual density-CMB two-point cross correlation measurement. We thus 
propose to measure the kSZ signal through the 0-0 cross correlation. This approach has a 
number of advantages. (1) Due to the underlying directional dependence of 0, it is uncorre- 
cted with the primary CMB, the thermal SZ effect and astrophysical contaminations such as 
the dusty star forming galaxies. Thus the 0-0 cross correlation picks up the kSZ signal in 
the SZ survey with a clean manner. (2) With the aid of galaxy redshifts, the cross correlation 
recovers the redshift information of the kSZ signal and allows for more detailed investigation 
on missing baryons. (3) Since the galaxy surveys usually have high S/N, the S/N of the kSZ 
measurement through the 0-0 cross correlation can be significantly improved. 

We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynam- 
ical simulations. We confirm that is indeed tightly correlated with at ^ < 1/z/Mpc, al- 
though nonlinearities in the density and velocity fields and nonlinear redshift distortion do 
weaken the tightness of the 0-0 correlation. We further quantify the reconstruction noise in 
from galaxy distribution shot noise. Based on these results, we quantify the applicability 
of the proposed kSZ tomography for future surveys. We find that, in combination with the 
BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to de- 
tect the kSZ with an overall significance of ~ 5Qcr and further measure its redshift distribution 
at many redshift bins over < z < 2. 

Key words: cosmology: observations - large-scale structure of Universe - cosmic microwave 
background 



1 INTRODUCTION 

Robust evidences from CMB and BBN show that the baryonic mat- 
ter accounts for ~ 4% of the total matter and energy of the universe. 
However, only a fraction of this baryon budget has been detected in 
the local universe, either in the form of stars, interstellar medium 
(ISM) and intracluster medium (ICM) (refer to lFukugita & PeeblesI 
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l2004l2006l for census of the baryon budget), while ~ 50% of the 
baryons remains elusive to robust direct detection. 

Looking for these "missing" baryons is crucial for the valid- 
ity of our standard cosmology model. The standard theory of hier- 
archical structure formation models indicates that the majority of 
baryons exist between galaxies as the diffuse intergalactic medium 
(IGM). Numerical simulations in the standard cosmology further 
suggest that a large portion of th is IGM is in the form of warm-hot 
intergalactic medium (WHI M) l lCen&Ostrike3ll999l :l Dave et al.l 
l200ltlCen & Ostrikedllooeh with temperature lO^K < T < 10^ K, 
which is believed to reside in moderately overdense structures such 
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as filaments. On one hand, this kind intergalactic medium is ionized 
dominantly by collisions, and is transparent to Lya radiation, and 
thus is hard to trace by Lycf forest. The absorption signature toward 
a bright X-ray source is also too weak to be resolved by current 
spectrographs. On the other hand, though WHIM emits radiation in 
the UV and soft X-ray bands, the emission strength is too weak for 
current instruments to detect. There have been observations trying 
to detect WHIM via absorption and radiation signature, but many 
of the m are of low significance or misinterpreted. See jBregmanl 
l2007h for a recent review of the status of WHIM detections. 

Among the methods of probing missing baryons, the Sunyaev 
Zel'dovich (SZ) effect is one of the most promising. Our universe 
is almost completely ionized after z = 6, where free electrons are 
prevailing in galaxy clusters as high energy ICM and in the less 
overdense filamentary structures as IGM. Free electrons will scatter 
off" CMB photons through inverse Compton scattering, and generate 
secondary CMB temperature anisotropics, which is known as the 
Sunyaev Zel'dovich effect. Therefore, the SZ effect is contributed 
by virtually all electrons. In principle, from the SZ observations, 
we are able to find all electrons and hence all baryons, due to the 
electric neutrality of the universe. 

The two major categories of the SZ effect are the thermal SZ 
effect (tSZ), arising through the thermal motion of electrons, and 
the kinetic SZ effect (kSZ), arising through the bulk motion of elec- 
trons. The efficiency for given electrons to generate the SZ effect is 
proportional to the thermal temperature for the thermal SZ effect 
and is proportional to the bulk peculiar velocity for the kinetic SZ 
efl^ect. Since the thermal temperature of electrons strongly couples 
to the electron density, the dominant contribution to the thermal 
SZ effect co mes from the high-d ensity and high-temperature ICM. 
For example. Iwhite et al.l ( l2002h found that 75% of the total ther- 
mal SZ effect at multipole { < 2000 comes from virialized regions 
with g as overdensity 6^^^ > 100 and lHernandez-Monteagudo et al.l 
( l2006l) found similar results: 80% of the tSZ signal comes from col- 
lapsed structures. Since the missing baryons are likely in less dense 
regions with lower temperature, the ability to find missing baryons 
through the thermal SZ effect is limited. 

On the other hand, the kinetic SZ efi"ect has a better poten- 
tial to probe the missing baryons. The peculiar velocity is deter- 
mined by the large scale gravitational potential and is thus only 
weakly coupled to local mass concentration. Hence, the contribu- 
tion to the kSZ effect is roughly proportional to the total mass of 
each baryon component, e.g., ICM and IGM. Namely, it is an ap- 
proximately unbiased probe of baryons, regardless of their thermal 
state (as long as they are ionized). Since the mass fraction of the 
missing baryons is ~ 50%, we thus expect a comparable contribu- 
tion to the kSZ effect from the missing baryons. This makes the 
kinetic SZ effect a promising probe of the missing baryons. For 
example, there have been works studying t he contribution of kSZ 
effect by WHIM to the CMB an isotropies jAtrio-Barandela et afl 
l2008l : lGenova-Santos et al.ll2009h . 

The remaining question is to measure the kinetic SZ effect, 
for which there are several obstacles. First of all, the signal of the 
kinetic SZ effect is not only weak, but also lack of spectral fea- 
tures to extract from the overwhelming primary CMB. Thus, from 
the CMB measurement alone, we only expect to detect it in the 
auto power spectrum measurement at f > 3000, where the pri- 
mary CMB da mps significantly and the kinetic SZ effect begins to 
dominate (e.g. lMa & Frvll2002l : Iziiang et alj|2004h . and at v = 217 
GHz frequency band, where the (non-relativistic) the rmal SZ effect 
vanis hes. However, refer to an interesting paper by iNozawa et al.l 
l2006l concerning the residual tSZ effect at 217GHz due to the rela- 



tivistic corrections by massive clus ters. Other conta minations such 
as the dusty star-forming galaxies jHall et al]|2010h make the kSZ 
detection even more difficult. For these reasons, even for CMB ex- 
periment as advanced as Planck, which has multiple bands over 
wide frequency range, it is difficult to detect the kinetic SZ effect 
directly, given its limited angular resolution. Secondly, like the ther- 
mal SZ effect, the kinetic SZ effect suffers from severe projection 
effect. It measures the electron momentum projected along the line 
of sight, thus the redshift information of baryons is entangled in 
the projection with a projection length of the size of the horizon. 
This not only degrades its power to infer the evolution of missing 
baryons over the cosmic epoch, but could also lead to large bias. 
For example, the patchy reioniza tion could contribute a significant 
fraction to the kinetic SZ effect l lZahnetalj|200^ : lMcO uinn et al.l 
l2005l : lmevet"ai]|2007h . Without redshift information, unless the 
reionization process is well understood and interpreted appropri- 
ately, the extra contribution from patchy reionization could be mis- 
interpreted as the sign of missing baryons. It's int eresting to notice 
the work of lHemandez-Monteagudo & Hoi | |2009|) . which proposed 
to recover the signature of the bulk flow of the missing baryons 
by cross correlating future CMB data sets with kSZ estimates in 
galaxy clusters. 

In this paper, we propose a kinetic SZ tomography method 
to overcome the above obstacles. The basic idea is to cross corre- 
late the CMB observation with a galaxy redshift survey or other 
surveys of the large scale structure with sufficiently accu rate red- 
shift information such as the 21cm intensity mapping tChang et al.l 
l2008h . Due to the cancellation mechanism arising from the under- 
lying directional dependence of the kinetic SZ signal as a result of 
its vector nature, direct cross correlation between the kinetic SZ 
effect and the galaxy number density vanishes (refer to Fig. [T] for 
more detailed explanation). This differs significantly from the tight 
correlation between the thermal SZ e ffect and the galaxy density 
jZhang & Penll200ll : lshao et al.l200^ . One way to circumvent this 
cancellation is to squa re the kinetic SZ effect and measure its cor- 
relation with galaxies jPore et alj2004lDeDeo et al.ll2005l) . 

We propose an alternative approach to avoid the cancellation. 
One can reconstruct the peculiar velocity field from the galaxy 
spectroscopic redshift survey, weigh it with the observed galaxy 
number density and other redshift-dependent factors to reconstruct 
a weighted momentum field. This weighted momentum field has 
roughly the same directional dependence as the true kSZ signal and 
thus avoids the cancellation. We thus expect a measurable cross cor- 
relation between the reconstructed field and the CMB map, which 
extracts the kSZ component in the CMB temperature fluctuation. 

The cross correlation signal should arise solely from the cor- 
responding redshift range where galaxies reside, to an excellent ap- 
proximation for a reasonably thick galaxy redshift bin (Az > 0.2) 
and sufficiently small angular scales (E > 10). Thus it disentangles 
the kinetic SZ contribution in this redshift range from contribu- 
tions of any other redshifts and recovers the redshift information 
of the kinetic SZ effect. It is for this reason that we dubbed this 
method as the kinetic SZ tomography, analogous to the thermal SZ 
tomography and also to the well known lensing tomography. An 
immediate application of the kinetic SZ tomography is to separate 
the patch reionization from the late time kinetic SZ effect. It is also 
effective to eliminate contaminations from the primary CMB, ther- 
mal SZ and other foreground contaminations, which do not have 
the characteristic directional dependence and thus should be uncor- 
related with the reconstructed momentum field. Since the galaxy 
surveys usually have excellent signal to noise (S/N), the cross cor- 
relation measurement can achieve much higher S/N than the auto 



The kinetic SZ tomography with spectroscopic redshift surveys 3 



CMB LSS 




Figure 1. The vanishing correlation between the kinetic SZ effect and the 
galaxy number density. The line of sight is from the observer (O), through 
the points B, E and A, to the last scattering surface. One galaxy resides at 
point G and we shall discuss the velocity distribution given the existence of 
this galaxy. The line EG is perpendicular to the hne of sight and AE=EB. 
The velocity at any point along the line of sight, e.g., the point A, can be 
decomposed, with respect to the position of the given galaxy, into one com- 
ponent along the direction AG and one perpendicular to the direction AG. 
The perpendicular component has equal probability to be along the direc- 
tion AC or the direction AD and thus the net contribution to the kinetic 
SZ-galaxy cross correlation is zero. So the only velocity component which 
may contribute is the one along the direction AG. However, from the sym- 
metry argument, if the evolution effect (light cone effect) can be neglected, 
velocity at point B has an equal probability to have a component along the 
direction BG, with the same amplitude. The projection of the two onto the 
line of sight cancels exactly. So the cross cotTelation between the kinetic SZ 
effect and the galaxy density vanishes. 



correlation measurement of thie kinetic SZ effect from tfie SZ sur- 
veys alone. Later we will show that, PLANCK plus BigBOSS can 
detect the kinetic SZ effect at ~ SOcr leveL 

At the beginning stage of this work. lHo et alj | |200^ published 
a work based on similar idea, whose applicability is further con- 
firmed in this paper. The two works are carried out independently 
and thus differ in many details. We test the proposed kSZ tomog- 
raphy against a controlled set of hydrodynamical simulations and 
quantify the tightness of correlations between the reconstructed 
map and the true kSZ signal at various redshifts over < z < 2. We 
further investigate several complexities such as redshift distortion 
and feedback, and demonstrate the robustness of the kSZ tomogra- 
phy against various complexities. Both works confirm the power of 
the kSZ tomography. 

The paper is organized as follows. In ^ we introduce the 
kSZ tomography method and discuss its limitations in general. We 
then test it against our hydrodynamic simulations in ^ The recon- 
struction and the test are done in both the real space 03.lt and the 
redshift space ( 33.21 ). In this section, we approximate galaxies as 
dark matter particles. So it corresponds to idealized surveys of vir- 
tually infinite galaxies such that shot noise in galaxy distribution 



is negligible. Realistic surveys have much lower galaxy number 
density and thus shot noise in the galaxy distribution induces non- 
negligible reconstruction noise. We directly quantify it from our 
simulations with a proper scaling 04.lt . We are then able to fore- 
cast its performance for survey combinations like the Planck CMB 
experiment plus the BigBOSS spectroscopic redshift survey (^. 
We discuss and conclude in ^ We present more technical details 
and further discussions in the two appendices. 



2 THE KINETIC SZ TOMOGRAPHY 

Bulk motions of free electrons induce secondary CMB 
anisotro pics, namely the kinetic S unyae v Zel'dovich (kSZ) 
effect l lSunvaev & ZeldovichI 1 1971 Il98(]|) . with temperature 
fluctuations 
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where 0, is the contribution from electrons in the i-th redshift 
bin, spanning the comoving coordinate range Xi ~ ^Xil^ < A" < 
Xi + A^,72. X is the comoving radial coordinate, c is the speed of 
light and t is the Thompson optical depth, p = (1 5s)v is the 
(normalized) electron momentum and the subscript "||" denotes the 
projection along the line of sight. 5^ is the electron number overden- 
sity, and v is the electron peculiar velocity. The weighting function 
Wksz = ^<.«eCrj- exp(-r)a/c modulates the contribution from each 
redshift. Throughout the paper, we focus on the kSZ effect after 
reionization and thus set^^ = I, an excellent approximation. Due 
to the neutrality, 5^ = (Jga,, where (5„as is the overdensity of (ionized) 
gas. 

As a vector, p can be always decomposed into a gradient (ir- 
rotational) part p^ and curl (rotational) part ps, p = p^ -I- ps, where 

V X p£ = and V • p^ = respectively. Here, the "E" and "B" 
notations are analogous to the electromagnetic fields. As can be in- 
ferred from Eq.[T] the contribution from the gradient part is largely 
canceled out when integrating along the line of sight, as long as the 
weighting function Wksz varies slowly across a correlation length 
of p£, which is of the order 100 Mpc//i todayQ. For the kinetic SZ 
effect after reionization, Wksz only changes significantly over the 
Hubble scale, thus the contribution from the gradient par t is negli- 
gible and the only significant contribution comes from ps jVishniad 
Il987l) . 

Pfi in general has two sources of contribution. Since V x pa = 

V X p = (1 -I- (5e)V X V + V(5p X V, the B-mode of p can come from 
the B-mode of v or from the cross talk between the density and 
velocity. Following the same notation, we can decompose the ve- 
locity into a "E" mode (gradient part) vg and a "B" mode (rotational 
part) Vg. For purely gravitational interaction, the velocity "B" mode 
decays, until multi-streaming and shell crossin g arise due to the 
nonlinear evolution (see for example chapter 2 of lBernardeau et ahl 
|2002i for a discussion). Thus in the linear and weakly nonlinear 
regimes a good approximation is v = and the only contribution 
to the kSZ effect comes from the cross talk between the density 



' This condition can be violated at the epoch of reionization, where the 
patchy reionization causes W to vary significantly over ~ 10/i/Mpc scales. 



4 Shao et al. 



gradient an d the velocity. This is the well known Ostri ker-Vishniac 
(OV) effect dOstriker & VishniacI igsdlVishniad 1987h . In the non- 
linear regime, Vg g rows and can also contribute to the kSZ effect 
dZhang et alj2004l) . 

The kinetic SZ tomography that we propose requires combi- 
nation of a SZ survey and a galaxy spectroscopic redshift survey 
with overlapping sky coverage. It contains three major steps: 

• Construct a 2D map 0; from the 3D distribution of galaxies 
in the i-th redshift bin. Ideally, 0, should be tightly correlated with 
0,, the true kinetic SZ signal from this redshift bin. A crucial in- 
gredient to guarantee a tight correlation is to estimate the peculiar 
velocity through the 3D galaxy distribution and use it to construct 
0,. Hereafter, we often neglect the subscript "i" where it does not 
cause confusion. 

• Cross correlate the reconstructed © with a overlapping SZ sur- 
vey. The cross correlation signal, to an excellent approximation, 
solely comes from the kinetic SZ in the chosen redshift bin. It is 
this step that recovers the redshift information of the kinetic SZ 
effect, eliminates various contaminations and reduces statistical er- 
rors. 

• Interpret the measured cross correlation signal and reconstruct 
the true kinetic SZ sign al. Tric k similar to the one adopted in the 
thermal SZ tomography l lShao e t al. 2009) can also be applied here. 

The current paper will focus on the first two steps and only briefly 
discuss the third step. 

2.1 The kinetic SZ reconstruction 

The primary goal of this paper is to extract the kSZ signal as well as 
its redshift information through cross correlating SZ surveys with 
galaxy spectroscopic redshift surveys. As explained early, the usual 
two-point cross correlation between the kinetic SZ effect and 
galaxy overdensity 5g vanishes {{®6g) = 0), due to the cancellation 
of positive and negative velocities along the line of sight (Fig. [T]l. 
A natural step to avoid such cancellation is to recover the velocity 
information and weigh the galaxies accordingly. This can be done 
in spectroscopic galaxy redshift surveys. 

Sectroscopic galaxy redshift surveys, such as 
LAMOSlH 
BigBOS^ SKAB Eucliciand JDEM/ADEPll], will mea- 
sure the 3D galaxy distribution ^j(x). At least part of them will 
overlap with SZ surveys such as AClB SPlH and PLANCK0 on 
the sky coverage. We are able to recover the velocity field from 
6g(x) and then 0, the galaxy momentum properly weighted. Since 
the direction of the peculiar velocity is taken into account in 0, the 
cross correlation (00) no longer suffers from the usual cancellation 
and thus <00> + 0. 

In the linear regime, the mass conservation (i.e. the continuity 
equation) reduces to 

<5,„ -H V • V = , (2) 
where 5„, is the matter overdensity and v is the peculiar velocity 

^ tittp://www.lamost.org/website/en 
^ http://cosmology.lbl.gov/BOSS/ 
http://bigboss.lbl.gov/index.html 
^ http://www.skatelescope.org/ 
* http://sci.esa.int/euclid 
' http://jdem.gsfc.nasa.gov/ 
^ http://www.physics.princeton.edu/act/index.html 
' http://pole.uchicago.edu/ 

http://www.rssd.esa. int/index.php?project=planck 



of the matter field. In the same regime, namely at large scale, the 
3D galaxy distribution is a good proxy of the underlying 3D mat- 
ter distribution. This can be described by a galaxy bias Sg = bgS,„. 
Under the above condition, given the observed 3D galaxy distribu- 
tion 6f'{x), we are able to obtain an estimator of the velocity field, 
which in Fourier space reads 

Kk) = -ifHSfHk)^ , (3) 

where It is the 3D wave vector, / = dlnD/d\na and D is the lin- 
ear density growth rate. The superscript "hat" in v (v) and in other 
symbols (e.g. 0) denotes the reconstructed quantity. In reality, to 
suppress the noise and stabilize the reconstruction, we often ap- 
ply some filters to the density field, before applying Eq. |3] So S"'"' 
should be treated as the smoothed density field. 

Eq.[3]is a biased estimator of the true peculiar velocity v. (1) 
The galaxy bias causes v to be overestimated by a factor bg. (2) 
Redshift distortion causes the observed density to deviate from the 
underlying matter density and thus biases the velocity reconstruc- 
tion. (3) Nonlinearities in the density evolution causes deviation 
from Eq. [2] (4) In the nonlinear regime where shell crossing and 
multi-streaming happen, velocity vorticity (rotational part, or curl 
part) develops, which is completely missed by the estimator Eq.[3] 

The imperfectness of the velocity reconstruction is not as se- 
vere as it looks, for the kinetic SZ tomography. Later we will show 
that the performance of the kSZ tomography is insensitive to deter- 
ministic errors in the velocity reconstruction. Thus a deterministic 
galaxy density and velocity bias, uncertainties in / and H, linear 
redshift distortion (the Kaiser effect) do not degrade the kSZ to- 
mography. However, stochastic errors from stochastic galaxy bias, 
nonlinearities in the evolution of density and velocity do. Later we 
will quantify their impacts and show them to be moderate at rele- 
vant scales. 

With the 3D density field Sg and thereby the recovered 3D 
velocity field v in hand, we can then reconstruct a weighted 2D 
momentum map 

& = J dxPnWiz), (4) 

where p = (1 + Sf'^)y{x), and v(x) is the inverse Fourier transform 
of v(k). We choose the weighting function W(z) = Wksz(z)- The in- 
tegral in Eq.|4]is over the corresponding redshift bin. Again, we can 
also apply some filters to the density field, before taking the prod- 
uct in p. These filters are not necessary to be the same as the ones 
for the velocity reconstruction. Furthermore, the density measure- 
ment in the product may not even be the same a s the one used i n 
the velocity measurement, as correctly pointed bv lHo etal.N2009t) . 

Here we want to clarify a likely confusing point. As discussed 
before, the kinetic SZ effect is mainly contributed by pg instead 
of p£. On the other hand, the reconstructed v in Eq. [3]is actually 
P£, which becomes clear later in Eq. [8] It thus seems that the re- 
construction misses the dominant contribution to the kinetic SZ ef- 
fect and thus should fail to work. However, there is an extra factor 
(1 + 6,) in the estiniator (Eq. [4). Recal l that, in the OV effect 
JOstriker & Vishniac|[l98^ : IVishniaclll987l) . it is the cross-talk be- 
tween the density gradient and the curl-free velocity that generates 
a curl component in p. Here, the cross-talk between the density 
and Pe generates a B-mode in the reconstructed p, which is tightly 
correlated with the true pg on relevant scales. This explains the 
reasonable performance of the reconstruction technique and the ki- 
netic SZ tomography. 
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2.2 How to quantify the kSZ tomography performance? 

The estimator is certainly imperfect. It can have both systemati- 
cal and random offsets with respect to 0. These deviations should 
be both scale and redshift dependent. These deviations can be vi- 
sualized by a ©-0 plot. Alternatively, it can be quantified by two 
parameters, r and h^. r describes the tightness of the 0-0 corre- 
latioiO and Z)g is the bias in with respect to 0. Clearly, both r 
and depend on the redshift range of galaxies used for recon- 
struction. We want to quantify the redshift dependence of r and b. 
Thus later in the analysis we will choose a projection length across 
which evolutions in r and be are negligible. For such projection, r 
and ig are functions of z and the 2D wave vector k±, the inverse of 
the perpendicular spatial separation. 



and 



Pe(k^,z) 



(5) 



(6) 



Later we will recognize as the perpendicular component of the 
usual 3D wave vector k. The Ps are the corresponding power spec- 
tra. 

There is no guarantee that b^ is close to unity, even in the lin- 
ear regime. As explained earlier, the velocity estimation is biased. 
Further, in the product oc (1 -i- S°^^)y, 6f'^ is also biased with re- 
spect to 6g. Later in the Appendix |X]we will show the complicated 
behavior of fog. 

However, large deviation of fog from unity does not necessar- 
ily mean poor performance of the kinetic SZ tomography, for three 
reasons. (1) First, the signal-to-noise ratio (S/N) of the (00) mea- 
surement, is solely determined by r. The S/N of each mode (Fourier 
or multipole mode) is 



P2. 
00 



^gg + (^0 + 0(^0+0 



1 + r-'- 



1 + 



(7) 



where Pg and are the corresponding noise power spectra. We 
find that the ig dependence drops out in the error estimation. A 
rescaling — » Zjg0 leaves no effect on the S/N of the cross corre- 
lation measurement, since both P'^,P@ fo? and r is unchanged 
under this scaling. (2) b^ t I definitely affects the expectation 
value of the cross correlation. However, since in the theoretical 
interpretation of the measured cross correlation, one can take the 
galaxy bias, redshift distortion and possibly other complexities into 
account and thus avoid systematical errors induced by fog 1 rea- 
sonably well . (3) Based on th e same technique in the thermal SZ 
tomography jShao et aL we can combine the cross correla- 

tion measurement Pqq and the auto correlation measurement Pg to 
obtain P@ = P~^l(r~P@)- This estimation relies on no information 
of fog and thus avoids the bias problem. 



^' r that we define differs from the one defined in lHo et'al First, 
their r is for the velocity field instead of the momentum field, 3D instead 
of 2D. Second, their r is not the cross correlation coefficient, but is actually 
analogous to rjb of the 3D velocity field, in our notation. 



For these reasons, in the main text we will focus on r to quan- 
tify the performance of the kinetic SZ tomography and leave dis- 
cussion on fog in the Appendix lAl 



2.3 Origins of the stochasticity r + 1 

A number of approximations made in the reconstruction pipeline 
cause the stochasticity in the 0-0 relation (r 5^ 1). (1) First of all, 
as an approximation to the exact mass conservation equation 



(5,„ + V • (1 + (5,„)v = , 



(8) 



the starting point Eq. [2] only holds where 5„, « 1 . (2) From 
Eq. |2] to Eq. [3] we have made assumptions of linear evolution 
{5,„{k,z) oc D(z)6i„(k, Zj)), deterministic bias between Sg and S,,,, 
and curl-free velocity. None of these approximations are exact in 
the nonlinear regime. (3) Even under these assumptions, Eq.|2]and 
[3] only hold in real space. Namely, we have implicitly assumed 
that the observed 6"^" is the galaxy overdensity Sg in real space 
(Sf'^ = Sg). However, in reality, what we directly measure is the 
galaxy number overdensity Sg in redshift space (S"^" = Sg). Due to 
the redshift distortion, Sg t Sg. The Sg-Sg relation is stochastic, due 
to nonlinear mapping between redshift and real space, nonlineari- 
ties in both the densit y and velocity field s (e.g. IWhite et"al]|2009l) 
and velocity vorticity jCarlson et al.l200^ . (4) Even if we have per- 
fect E-mode velocity reconstruction, we still have no handle on the 
B-mode velocity, which contributes the kinetic SZ effect. (5) In Eq. 
[4] we multiply the reconstructed velocity with 1 -1- 5°'"' instead of 
I + Se. The possible stochasticity between Se and Sg also increases 
the stochasticity between and 0. 

With the aid of our hydrodynamic simulation, we are able to 
quantify the combined influen ce of all t hese f actors on r, except 
for the stochastic galaxy bias. IHo et alj | |2009|) used the halo oc- 
cupation model and N-body simulations to produce galaxy mock 
catalog. This approach captures the galaxy stochasticity and shows 
that the kSZ tomography is robust against it. Our simulations have 
relatively small box size (lOOfo^'Mpc) and hence do not allow us to 
follow the same approach. So we leave this issue for further inves- 
tigation. 



3 TESTING AGAINST HYDRODYNAMICAL 
SIMULATIONS 

We test the kinetic SZ tomography against our hydrodynamical 
simulations. Th e simulations are run with the GADGET2 code 
jSpringell |2005|) in a ACDM cosmology with parameters: A = 
0.732, no = 0.268, n,, = 0.044 h = 0.71, o-g = 0.85. The box 
size of the simulation is L = lOO/i"' Mpc on each side, in which 
512' dark matter particles and 512' gas p articles are initially seeded 
(See more details about the simulation in Ijing et al.l2006l ; lLin et al.] 
I2OO6I) . We have an non-adiabatic run, in which gas particles are al- 
lowed to cool and condense into collisionless star particles, along 
with which SN feedback is taken into account. We also have an 
adiabatic run with the same cosmological parameters and the same 
initial conditions. We will focus on the non-adiabatic run, since it 
has better capture on the gastrophysics and hence better modeling 
of the kinetic SZ effect. Unless otherwise specified, all simulations 
are based on this simulation. We also analyze the adiabatic run to 
better understanding the generality of the kSZ tomography ( 33.3.2t . 
We choose in this work a few representative redshifts to quantify 
the feasibility of the tomography technique, which is mainly char- 
acterized by the quantity r(k^,z). 
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Our hydrodynamic simulations hiave direct information of the 
gas momentum distribution and thus the kinetic SZ effect. How- 
ever, the simulations don't simulate galaxies. To proceed, we ap- 
proximate galaxies as simulation dark matte r particles. Since the 
galaxy stochasticity is likely sub-dominant ( lBonoH&Pe"i]|2009l ; 
iBaldauf etal.l201(i) . this approximation is reasonable to estimate r 
between and Q, althoug h it indeed over-estimates it. As demon- 
strated bv lHo et al] | |2009|) . the kSZ tomography is robust against 
the galaxy stochasticity. So we will leave this issue elsewhere. 
Since the number density of simulation particles is much higher 
than that of galaxies in any realistic surveys, shot noise in the 
reconstruction is negligible. To use the measured r for forecasting, 
we have to take the shot noise inevitable in realistic surveys into 
account. We leave this issue until next section. 

We test the kinetic SZ tomography at several typical redshifts 
z = 0,0.53,1.02,2.08. We reconstruct the velocity and conse- 
quently the momentum field from the dark matter distribution in the 
corresponding simulation outputs. At these redshifts, we project the 
3D momentum field over a single box size (100/;"' Mpc ) to get the 
2D momentum maps. We then compare these maps with the maps 
of underlying kinetic SZ signal, which is directly measured by mul- 
tiplying the gas density and velocity distribution in the correspond- 
ing simulation outputs. Since the simulations have information of 
peculiar velocity, we carry out the reconstruction and comparison 
in both real space and redshift space. 



3.1 The reconstruction in real space 

3.1.1 The reconstruction procedure 

We first carry out the reconstruction in real space. The reconstruc- 
tion is a two-step process. The first step is to reconstruct the pe- 
culiar velocity. At each redshift z,, we construct the matter density 
field (5(x) using clouds-in-cells (CIC) scheme in a partitioning of 
= 256^ cells, followed by a Fourier transform to get the matter 
distribution in Fourier space 5(k). We then apply a Gaussian fil- 
te{3 ^a(k) = sxp{-k^R^J2) to the original field i5(k) in order to 
wipe off nonlinear fluctuations on small scales, which are uncorre- 
lated with the large scale velocity. The adopted smoothing length 
Rs = 1.56/;"' Mpc applied in the whole context is around the radius 
of a typical cluster, and we will discuss the influence of different 
smoothing length in 33.3.11 We then obtain the reconstructed ve- 
locity v(k) from Eq.[3]through this smoothed density field. By the 
inverse Fourier transform, we obtain the real space velocity field 
v(x). 

The second step is to project the momentum field along the 
line of sight h to get the reconstructed 2D map at redshift z, 



0;ij2(n) = 



(9) 



Since the weighting function Wksz changes little over that scales, 
so we treat it as a constant and omit it. The sum is over a single 
simulated box at the investigated redshift, instead of stacking the 
box in the light cone. Notice that the density term in the above 



IHo et al.l j2009l) populated halos with galaxies. Since in that case the 
galaxy density is low, and shot noise is non-negligible, they applied the 
Wiener filter to reduce the Poisson noise. This is also a necessary procedure 
when dealing with real data. In our work, we use simulation particles as 
galaxies, whose number density is much higher and thus the shot noise is 
negligible. For this reason, we instead apply generally a Gaussian window 
function to filter out small scale nonlinear fluctuations. 




Figure 2. The 0-0 map (in arbitrary unit) in the real space at z=0, 0.53, 
1.02, 2.08. In order to further reduce the scattering of 0-0 relation, we av- 
erage the recovered velocity field in every neighboring 4^ = 64 cells to get 
the averaged velocity field v as well as the density field i5 on a much coarser 
partitioning, say, 64^ cells. We show here the cell-to-cell correspondence in 
all three Cartesian directions by randomly selecting 512 out of 64^ couples. 
The recovered field intimately follows the kSZ signal along the diago- 
nal at high redshifts, while at low redshift, there are great fluctuations due 
to non-linear evolution of the density field and thereby the velocity field. 
The three components in Cartesian directions are denoted by different point 
types, triangles for x, squares for y and circles for z. 



equation (Eq. [9} is the unsmoothed density, which preserves the 
necessary information of the true density field in the deriving of 
the momentum field. In practice, we construct at each redshift three 
maps along three Cartesian coordinate directions of the simulated 
box respectively. As they can be considered as three independent 
measurements, in the following figures which concern statistics, we 
show the average results unless specified. 

At the same time, we obtain the true kSZ signal in a similar 

way 



0;ij2(n) = 2 n 



. '^bji.p-.fl 



(1 + Sbjljl.fi) : 



(10) 



where Si, is the nonlinear overdensity of baryons. Vj, is the bulk 
velocity derived directly from the simulation by averaging within 
the host cefl (jl,j2,j3) 



(11) 



where v, and nij are the comoving peculiar velocity and the mass of 
i-th gas particle in the host cell, w, is a spline kernel used to smooth 
the gas particles. 

3.1.2 The @-@ relation 

In Figure|2]we show the cell-to-cell 0-0 correspondence. The 0-0 
data points scatter around = 0, meaning a bias bg ~ I. However, 
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Figure 3. The cross correlation coefficients r between and 0. The cor- 
relation coefficient is averaged over the three components. The two are 
highly correlated, r=:0.9, on lai'ge scale kx_ S5 1/iMpc^'for all the 4 red- 
shifts, while less correlated to smaller scales ends due to the non-linearities 
coming into play. Nevertheless, the correlation coefficient is still tight at 
z=0, with r > 0.6 up to kj^ ~ Hit Mpc"' . 



this should not be over-emphasized since this only represents an 
surrealistic case of galaxy bias bg = 1 without redshift distortion. 
Reconstructions based on galaxy sample with bg t 1 would result 
inbe i= I and hence a different slope of the 0-0 relation. Redshift 
distortion also changes the slope of the ©'-0 relation, as can be 
seen from Fig.|4] 

The 0-0 relation shows non-negligible dispersion around the 
mean, but still reasonably tight. This means that the stochasticity is 
noticeable, but not yet overwhelming. We also notice that, the dis- 
persion gets stronger at lower redshift. This is not surprising, since 
the reconstruction is based on linear theory and thus works better at 
higher redshift. Nonlinearities in the density field and velocity field 
degrades the reconstruction accuracy, as discussed in 32.31 



3.1.3 The cross correlation coefficient r 

Out of the two quantities concerning to the reconstruction perfor- 
mance, the cross correlation coelficient r describes the tightness of 
0-0 relation and is a major measure of the kinetic SZ tomography. 

is of less importance in quantifying the kSZ tomography per- 
formance, so we leave related results to the Appendix [A] To mea- 
sure them, we perform 2D Fourier transforms of (0). 0(kx) = 
Jt/-Xj^0(xj^)exp(ikx ■ Xj^)/A, where A is the area of the map. Note 
here and Xj^ are both 2D variables. We then obtain the power 
spectrum P@ with (iTrfdoikj. - kj_')Pe{kJ = <0(kx)0(kx')>- 

Fig. [3] shows that the stochasticity in the 0-0 relation is in 
general not a severe issue, even in the strongly nonlinear regime, 
consistent with Fig. [2] As shown we have strong correlations, r > 
0.9 to = 1/7 Mpc"' at all redshifts, and even r > 0.8 to kj^ ~ 
3/! Mpc" 'except z=0. Nonlinearities do degrade the reconstruction, 
as we see that r decreases towards low redshifts. However, even 



Figure 4. The 0-0' relation in the redshift space. The same routine as Fig. 
[2]is used in deriving this relation. At high redshifts, the recovered 0' over- 
weighs the underlying kSZ signal 0, because the Kaiser effect enhances 
the density by a factor ~ 1 + //3, and the higher redshift the higher the 
slope is, since / ^ Djjj^d) is larger at early epochs. On the other hand, the 
scatters in 0-0' is much larger than that in the real space due to the larger 
stochasticities of density-velocity relation and the real space-redshift space 
mapping. 



at z = 0, the 0-0 correlation is still pretty tight, with r > 0.6 to 
k± ~ SAMpc"'. We don't correct for the aliasin g effect, such that 
small scale results may be misleading ( |jind2005l) . Nevertheless, we 
do expect good results up to a quarter of the Nyquist wavenumber, 
i.e. around 2h Mpc"' , and it's safe for us to estimate the cross power 
spectrum up to f ~ 2000 - 3000. 



3.2 The reconstruction in redshift space 

In reality, what we get from a galaxy survey is the galaxy num- 
ber density in the redshift space. As expected in ^ the redshift 
space distortion will induce new sources of uncertainty and fur- 
ther degrades the reconstruction. We are able to quantify its impact 
through our simulations. Due to its own peculiar motion, the appar- 
ent position of each dark matter particle along the line of sight can 
be written as 



: X -I- ■ 



H 



(12) 



where n denotes the unit vector along the line of sight, and v is 
the comoving peculiar velocity, x' is position in the redshift space, 
and from here on the superscript ""' indicates the corresponding 
quantity in the redshift space. As a result of this displacement, what 
we observed is a distorted density field S". With as the starting 
point, we can reconstruct the momentum fields 0' through Eq. [3] 
following the same procedure as in the real space. 
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Figure 5. The 0-0' relation in tlie redsliift space after coirecting for tlie 
Kaiser effect. It's interesting to find that the slopes approach unity at all 
investigated redshifts, while the scatters keep nearly unchanged compared 
to Fig.E] 



3.2.1 0-0' relation in the redshift space 

The 0-0' relation is shown in Fig.|4l which is significantly different 
from the one in real space (Fig. [J}. We can see that in two aspects. 
(1) The average slope of the 0-0 relation changes, from to 
fl0 with fl > 1. This is caused by the linear redshift distortion 
(the Kaiser effect), which induces 

d'(^) = 5m\+PiilX (13) 

where ji^ = n • k is the cosine of the angle between the line of 
sight n and the wavevector k. j8 = f/bg. The Kaiser effect enhances 
the galaxy overdensity by a factor 1 -I- j8/3 > 1 and thus causes 
the slope a > I. We also notices that the slope a decrease from 
high redshift s to low redsh ifts. There are two causes. First, p = 
/ ^ n|;,«(a) jPeebles| [T980h (in our case where bg = I) decreases 
with decreasing redshift. Second, the finger of God effect caused 
by small scale random motion suppresses the redshift space galaxy 
density. This effect becomes stronger at lower redshift. (2) Scatters 
in the 0-0' relation are significantly larger than that in the real 
space (Fig.|2j. Stochasticities in the density-velocity relation and 
real space-redshift space mapping are largely responsible for these 
larger scatters. 

It's interesting to show the 0-0' relation if we correct for the 
Kaiser effect. The figure is shown in Fig.|5] We find that the slopes 
approach unity at all investigated redshifts, while the scatters keep 
nearly unchanged compared to Fig. |4] This is, however, not un- 
expected, since the deterministic Kaiser formula only changes the 
amplitude and is reversible given cosmology. Thus finger of God 
should be mainly responsible for the scatters of the relation. 



Figure 6. The cross coirelation coefficients r(kj^,zf in the redshift space. 
Thick lines show r' with Kaiser effect corrected for. Compared to Fig. [3] 
i'(f<^±,zy is suppressed throughout all scales, especially at k±^ > 1/jMpc"'. 
At I < 1, the reconstruction works poorly, with r ^ 0.3 on scales k±^ > 
3h Mpc"' . However, there's a considerable correlation strength r' > 0.5 for 
k± < IftMpc^'at z=0. At higher redshift where nonlinearity of velocity is 
not significant, the correlation strength is still tight enough, e.g. with r 2: 0.7 
forA-^ ~ l/!Mpc"'atz=2.08. 



3.2.2 r'(k^, z) in the redshift space 

The cross correlation coefficient r in redshift space is shown in 
Fig. |6] Compared to Fig. |3] the correlation coefficient is sup- 
pressed throughout all scales, especially at k^^ > l/7Mpc"'. At 
kx^ > 3/iMpc"'and z < 1, the kinetic SZ reconstruction and thus 
the kinetic SZ tomography works poorly since the reconstructed 
barely resembles the true signal @ (r < 0.3). These results show 
unambiguously that redshift distortion is a significant source of er- 
ror in the kinetic SZ tomography. It's also worthwhile to see the 
changes if we correct for the Kaiser effect. They're also shown 
in Fig.|6]as thick lines. As indicated by Fig.|5] correcting for the 
Kaiser effect does not influence r" too much on scales of interest. 

As explained in 32.31 the degradation is caused by the stochas- 
ticity in the 6-6" relation. They can be induced by the nonlinear 
mapping between the real sp ace density and the redshift space den- 
sity (e.g. lScoccimarrol2004l) and the stochastic ity between the den- 
sity and velocity field (e.g. IWhite et"al]|2009l) . This stochasticity 
then induces the stochasticity in the reconstructed velocity v, with 
respect to true velocity (Eq.[3]l. In the momentum reconstruction we 
need to multiply the reconstructed velocity by an extra factor 1+6 
(Eq. |4j to obtain the reconstructed momentum. So the stochastic- 
ity in the momentum reconstruction and thus the kinetic SZ recon- 
struction has an extra source, from the term \ + 6. All these com- 
plexities worsen the reconstruction, increase scatters in 0-0 and 
decrease r. 

Despite the above degradations, the reconstructed still 
shows reasonably tight correlation with the true signal at < 
1/iMpc"'. We will show in the next section that this correlation 
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Figure 7. The deviations between different smootliing scales, say, 0, 
0.78/r' Mpc, 1.56/1"' Mpc and 3.12/1"' Mpc , in both real space (top panel) 
and redshift space (bottom panel). In the real space, smoothing length with 
0.78/r' Mpc works best at z=2.08 while 3.12/;"' Mpc is the most suitable 
at z=0. However, in the redshift space, the enhancement due to smooth- 
ing is less effective, and a larger smoothing length is better. Probably an 
anisotropic smoothing other than a Gaussian smoothing should work better. 
As a median smoothing length, 1.56/!"' Mpc is effective at both low and 
high redshift. 



strength allows for robust kinetic SZ tomography, combining the 
Planck CMB experiment and the BigBOSS galaxy spectroscopic 
redshift survey, or other surveys with comparable power. 

An interesting behavior to notice is that at z = 0.53 is worse 
than at z = 0. This result is consistent with a larger scatters in the 
0-© plot at z = 0.53 than at z = 0, implying that stochasticities 
induced by the Finger of God effect becomes the largest at that 
epoch. At redshift higher than z ~ 1, the reconstruction suffers less 
from the nonlinearities and is thus stronger, e.g. with r > 0.7 at 
k± < 1/i/Mpc at z=2.08. Actually, we find in simulations that the 
comoving velocity dispersion peaks at z ~ 0.6. This would induce 
the largest anisotropics and hence make r' at z ~ 0.6 the worst. 



3.3 Uncertainties in r 

In the above section, we quantify the impact of redshift distortion 
on r and thus show that redshift distortion degrades the kinetic 
SZ reconstruction significantly. There are other factors affecting 
the reconstruction. Here, we briefly discuss the influence of fil- 
ters adopted to smooth the density field 03.3.U and gastrophysics 
(EH). 



3.3.1 The influence of density smoothing on r 

The results shown above all adopt a Gaussian filter with smoothing 
length = 1.56/r' Mpc to smooth the density field before velocity 
reconstruction. The reconstruction robustness no wonder depends 
on the way to smooth the density field. We do not aim to perform a 



k/h Mpc" 



Figure 8. The difference of the coirelation coefficient between adiabatic and 
non-adiabatic simulations. The overall r' is insensitive to gastrophysics, and 
there're only ~2-3% changes of r\ It will accordingly influence the estima- 
tor of the signal-to-noise of the cross power spectrum by several percent 
level, however, for kSZ tomography, this is sufficient for the current sur- 
veys. 



comprehensive investigation on this issue. Rather, we will restrict 
to the Gaussian filter and investigate the dependence of r on R,. We 
arbitrarily compare between the cases of = 0,0.78, 1.56,3.12 
/;"' Mpc . For clarity, we only show the comparisons at z = and 
z = 2.08 inFig.|7]. 

The basic (and obvious) conclusion is that, smoothing is nec- 
essary to suppress small scale nonlinearities and improve the re- 
construction. Fig. |7] (upper panel) shows general improvement in 
r when smoothing is taken, comparing to the case of no smooth- 
ing (Rs = 0). For example, smoothing the density field with 
Rs = 1.56/j"' Mpc can boost r by ~ 30% at k = l/iMpc"'and a 
factor of 2 or more at smaller scales, in real space. 

Gaussian smoothing is less effective in redshift space: the im- 
provement in Tj is often less than 10% on scales of interest (bottom 
panel, Fig.|7]l. The redshift space overdensity is anisotropic, so a 
spherical Gaussian smoothing will not work well. An anisotropic 
filter may work better in redshift space. This is certainly an inter- 
esting technical issue for further investigation. 

A right smoothing should balance between suppressing small 
scale nonlinearities and preserving large scale signal. If R, is too 
large, it may wipe off too much large scale clustering responsible 
for peculiar velocity and thus degrade the reconstruction (decrease 
r). This may be the reason that R^ = 0.78/j"' Mpc works better than 
larger R„ at z = 2.08 in the real space. However, overall R^ = 1.56 
/?"' Mpc works reasonably well at all redshifts investigated. 

The simulation data we deal with is like an ideal survey, with 
negligible shot noise, uniform selection function and regular survey 
boundary. Smoothing for real data is of course much more compli- 
cated. For example, a big issue in real survey is shot noise due to 
low galaxy number density, especially in spectroscopic redshift sur- 
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veys. For this issue, one can refer to lHo et alj | |2009|) for discussion 
on the application of the Wiener filter. 



3.3.2 r and gastrophysics 

r also depends on gastrophysics. All the results shown above are 
based on our non-adiabatic simulation, with radiative cooling, star 
formation and supernova feedback. Although we are not able to 
robustly quantify its detailed dependence on these gastrophysical 
processes, we can obtain a rough estimation by comparing the 
above results to our adiabatic simulation with identical initial con- 
ditions. The relative differences of r' between the two simulations 
are shown in Fig.[8] We find that, overall r and thus the performance 
of the kinetic SZ tomography is insensitive to the gastrophysics. r 
at fcx < 1/zMpc^'only varies by less than 2-3%. The influence of 
gastrophysics is larger at smaller scales, but is still less than 5% up 
toitx = 3/!Mpc-'. 

This insensitivity to gastrophysics has two implications on our 
kinetic SZ tomography. First, it is unlikely that some realistic gas- 
trophysical process not included in our non-adiabatic simulation 
can dramatically suppress r and thus invalidates the tomography. 
Second, this significantly simplifies the theoretical inter pretation of 
the to mography results. Based on the same technique in jShao et al] 
l2009l) , we can combine the two measured correlations, (00) and 
<00>, to obtain (00) arising from the same redshift bin. The only 
unknown quantity in this approach is r, which we shall calibrate 
against simulations. However, if r is very sensitive to gastrophysics, 
the calibration on r would be very difficult due to large uncertain- 
ties in our theoretical and numerical understanding of these gas- 
trophysics. The insensitivity of r to gastrophysics shown in Fig. 
[8] implies that, despite imperfect theoretical and numerical under- 
standing of these gastrophysics, r can still be accurate to a few per- 
cent level at relevant scales. This precision suffices for the kinetic 
SZ tomography. 



4 ERROR FORECAST 

The r measured above quantifies the performance of the kinetic SZ 
tomography for a virtually ideal galaxy survey, for which the simu- 
lation particle number density is high and the shot noise is negligi- 
ble. However, in real surveys, the galaxy number density is a factor 
of 10"* -10^ smaller, resulting in much larger and thus non-negligible 
shot noise. The shot noise affects both the velocity and momentum 
reconstruction. We use our simulation to quantify these effects and 
then apply these results to forecast the performance of the kinetic 
SZ tomography. Our target CMB experiment is Planck and the tar - 
get galaxy spectroscopic survey is BigBOSS dSchlegel et al.l2009h . 
The kinetic SZ tomography based on other SZ surveys like SPT and 
galaxy surveys like ADEPT, Euclid and SKA is expected to work 
better. 



4.1 Estimating the reconstruction noise from galaxy 
distribution shot noise 

The discreteness of dark matter particles induces spurious fluctua- 
tions (shot noise) in the density field, which we denote as 6^, and 
thereby in the reconstructed velocity, which we denote as Vjv. The 
correlation of the contaminated momentum 0^. is (for simplicity. 
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Figure 9. The ratio P^^ , /^0i °f 'he shot noise power spectrum and the 
reconstructed momentum power spectrum. Both the A^^ and A^^ are much 
smaller than the reconstructed momentum powers, with Ni terms domi- 
nating A^i towards smaller scales. Given a small number density in cur- 
rent galaxy redshift survey, the total shot noise would be comparable to 
the signal. Nevertheless, on scales of interest for kSZ tomography, say 
< 1/iMpc-', galaxy surveys would provide a relatively less contami- 
nated reconstruction of the momentum. 



we omit the projection along the line of sight) 

= {(i+6')r ■h(i+s")r' -h'} 

+{SnV ■ hS'^V ■ h') . 

The last two terms in the above equation are the leading non- 
vanishing noise terms. We denote the first noise term as oc 
(1 + (5')vn while the second as Ni °^ Sn'^' 

By randomizing the particle positions in the simulation, we 
can directly measure Sn and derive therein \n oc 6N{k)k/k. Fol- 
lowing the same procedure as in reconstructing the momentum 
map, we produce two corresponding maps of these two shot noises. 
In order to estimate the contribution of shot noise, we compute 
the power spectrum of the two noise terms, and show the ratios 
Pn, IP&s in Fig. |9] We can see that both the Pf^^ and Pjv, noise 
terms are much smaller than the reconstructed momentum powers. 
Although A'l and A'2 terms are comparable on the very large scales, 
they deviate significantly towards smaller scales, with A'2 terms pre- 
dominantly overweighing N\ terms. On the small scales, A'2 terms 
would turn out to be several thousandth of the reconstructed mo- 
mentum 0^, and they can be comparable to or even dominate the 
signal, given a much smaller number density in current galaxy red- 
shift surveys. Nevertheless, on scales of interest for kSZ tomogra- 
phy, say < 1/7 Mpc-', galaxy surveys would provide a relatively 
less contaminated reconstruction of the momentum. 
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4.2 Error forecast in C/ 

Based on the above results, we are able to estimate the S/N for real 
surveys. Since in this section we consider the error forecast in real 
survey, i.e. in the redshift space, we omit the superscript for 
simplicity. For a given redshift bin e [z - Az/2,z + Az/2], we are 
able to measure the angular cross power spectrum Cf between the 
reconstructed and the CMB measurement where the kSZ effect 
is embedded. The statistical error in this measurement can be esti- 
mated by 



Ci 



fCMB , f-kSZ , rCMB,N (-N 

1 + h. ' ' n + -i-) 

i ^ 9^kSZ.Az X'^ ^ c- ' 



(14) 
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Here, is the angular power spectrum of primary CMB and 

^CMB.N jjjg power spectrum of the measurement noise. The ther- 
mal SZ effect is also a source of noise. Since we work with the 
frequency band around 217 Ghz, where the thermal SZ effect virtu- 
ally vanishes, we will neglect the thermal SZ effect. There are other 
possible sources of error, such as the dusty star forming galaxies 
and radio sources. They are not likely dominant over the primary 
CMB at f < 2000 of our interest. Thus we will not include them in 
the error analysis. Cf''^ is the angular power spectrum of the kinetic 
SZ effect and C'f''^''^ is the contribution from the given redshift bin. 

/ Cq is the ratio of the angular power spectrum of the reconstruc- 
tion shot noise and of the reconstructed momentum 0. r is again the 
cross congelation coefficient between the reconstructed kSZ map 
and the true kSZ signal. However, now the projection length is in 
general much larger than the simulation box size 100/;"' Mpc. For 
interesting angular scales of ^ ~ 10^, two redshift bins separated 
by more than 100/7"'Mpc can be approximated as uncorrelated. So 
we can average over r measured in the last section over the relevant 
redshift range to obtain an estimation of the r used in this section. 

The first factor 1 in the rh.s of Eq. [14] comes from the cos- 
mic variance in the cross-correlation signal and assumes the ki- 
netic SZ effect and the reconstructed are Gaussian. The ac- 
tual non-Gaussianties will increase the cosmic variance. However, 
this cosmic variance term is sub-dominant to the other term, since 
(C™'' -nCf 2-HC™^ '')/Cf ^-^^ » 1 and r s; 1. So the Gaussian as- 
sumption is reasonable. And the second relation holds since usually 

(C™s cf^ + c™''-'')/cf ^-^^ » 1. 

From this equation, we can figure out the improvement in S/N 
of the kinetic SZ measurement. To better see this point, we will 
discuss under the limit that the galaxy number density is sufficiently 
high such that C^/Cg, <c 1 at sufficiently large angular scales. We 
then have 
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(15) 



We can compare it to the S/N of the kSZ auto correlation power 
spectrum measurement 



.CMB.N 



(16) 



As Cf^^ » Cf 2 at f < 2000, the improvement in S/N by the kSZ 
tomography is of the order r ^Cf^/Cf^ yjcf^-^^Cf^. Since r 
is close to unity at relevant scale (e.g. £ ~ 10^ and z ~ 1), for 
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Figure 10. The error of the measurements of cross power spectrum between 
at redshift bins z=0.2 - 0.6, 0.6 - 1.0 and 1.0 - 2.0, using BigBOSS plus 
PLANCK. There are significant signals around f ~ 2000, with en'ors 26%, 
li% and 6% respectivety for the thi'ee redshift bins. With a larger redshift 
bin and better r' over the reievant muftipofe ranges, the totaf S/N in z = 
1.0 - 2.0 is 46, the more considerabfe than the other two, fO at z=0.2-0.6 
and 21 atz=0.6-1.0. 



2 or beyond), the 



'™B/cf^ » 1. In reality. 



a sufficiently deep galaxy survey (e.g. to z 
overall improvement is of the order ^jc^ 
since the galaxy number density is finite, shot noise dominates at 
small scales, as can be inferred from Fig.|9] 

Now we numerically calculate the S/N for the combination 
of PLANCK plus BigBOSS-N. Cf^^ is calculated from CAME0 
The instrument noise power spectrum of PLANCK is 



CMB,N 



(17) 



exp(-f(^ + \)I2CI,J for a 
For PLANCK, 



C 

where the window function W[ 
Gaussian beam with ^beam = V8 In2/6'FWHML 
^FWHM = 5.0 arcminutes. Within which the average Icr sensitiv- 
ity per pixel CTp j/TcMB = 4.8 x 10"* is expected after 2 full sky 
surveys. PLANCK will survey the sky at several frequency bands. 
The 217 GHz band where thermal SZ signal vanishes is the one 
that we use to cross correl ate with galaxies) Cf ^ and cf ^■'^^ are 
calculated by the model of lZhang et al.l 120041) . 

For the target galaxy s pectroscopic redshift survey, we choose 
the BigBOSS(-N) project Ischlegel et al]|2009l) . It plans to mea- 
sure the spectroscopic redshifts of 30 million LRGs and emission 
line galaxies at 0.2 ^ z ^ 2 over 14000 deg^ sky coverage. Since 



http://camb.info/ 

http://www.rssd.esa.int/index. php?project=planck 
The 217 GHz band is the only band promising to detect the kinetic SZ 
effect in auto correlations. However, the kSZ tomography is able to detect 
the kSZ effect in other frequency bands and the resuiting S/N can be com- 
parable at= C < 2000 where the thermai SZ effect is sub-dominant to the 
primary CMB. This is an issue for further investigation. 
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PLANCK will survey the full sky, the overlapping fractional sky 
coverage is /^i^y = 14000/(4;r/(7r/180)^). Given the much lower 
galaxy number density (comparing to simulations), shot noise is a 
significant source of error in the kinetic SZ tomography. The result- 
ing reconstruction noise can be roughly estimated through a simple 
scaling 



where, Ps are the ones at some intermediate redshift of the redshift 
bin and k± = CJx- '^Am and Sobs are the surface densities in the 
simulation and in the real survey respectively. This estimation is by 
no means exact. But it is reasonably good to demonstrate the power 
of PLANCK plus BigBOSS to detect kSZ. 

We divide the BigBOSS survey into three redshift bins, 0.2 ^ 
z ^ 0.6, 0.6 ^ z ^ 1.0 and 1.0 ^ z ^ 2.0. This choice is somewhat 
arbitrary. In the error estimation, we simply adopt r and / Cg, es- 
timated at some intermediate redshift of each redshift bin. In the 
error estimation using Eq. [141 the multipoles { are logarithmically 
equal spaced, say, Af/f = 12.2%. We present the error distribu- 
tion of the cross power spectrum in Fig.[TOl We find that the cross 
correlation can be measured robustly at ^ ~ 1000-3000. For exam- 
ple, the relative errors at f = 2000 are about 26%, 11% and 6% at 
the three redshift bins correspondingly. Good performance at these 
scales is expected for three reasons. (1)0 and the underlying are 
tightly correlated at corresponding scale fc^ = ^Ix- (2) Shot noise is 
low for a galaxy survey as big as BigBOSS, comparing to the noise 
level in the kSZ measurement. (3) The primary CMB drops while 
the measurement shot noise is still bearable at ^ = 2000. 

The accuracy of the cross power spectrum measurement in- 
creases with redshift. This is partly due to stronger r at higher 
redshift. And this epoch, the comoving scale = (/x(z) moves 
to more linear regime and the nonlinearity becomes weaker. Both 
push r upward and improve the kSZ tomography. What's more, the 
binned kSZ signal from z= 1.0-2.0 is larger than the other two on 
scales of interest. It certainly helps the cross correlation measure- 
ment. 

On the other hand, it is unlikely to detecting kSZ at ^ < 300 
and C > 3000 through the kSZ tomography with PLANCK plus 
BigBOSS-N . Toward larger scales, the kSZ signal decrease, while 
the primary CMB overweighs by 3 or 4 orders. On smaller scales, 
the finite angular resolution of the PLANCK survey causes the 
CMB measurement noise to increase exponentially. At the same 
time, the shot noise in BigBOSS begins to dominate over the galaxy 
clustering signal. Here we want to caution the readers that the S/N 
at ^ > 3000 is overestimat ed, since w'e neg lected noises from the 
dusty star forming galaxies jHall et alj2010l) . which may dominate 
over the primary CMB. 

We estimate that, the overall S/Ns for the z = 0.2 - 0.6, 0.6- 1 .0 
and 1.0-2.0 redshift bins are 10, 21 and 46, respectively. The com- 
bined S/N of the whole redshift range z = 0.2 - 2.0 is 51. These 
results are consistent with the findings of lHo et ^ | |2009|) . in which 
they use PLANCK plus SDSS and ADEPT. Although these num- 
bers are likely only accurate within a factor of 2, they nevertheless 
robustly confirm the applicability of the proposed kSZ tomogra- 
phy to realistic surveys. Since the overall S/N~ 50, we are able to 
choose finer redshift bins and thus measure the kSZ evolution over 
~ 10 redshift bins. 



5 DISCUSSIONS AND CONCLUSIONS 

In this paper, we proposed the kSZ tomography method, tested it 
against simulations and estimated its applicability to realistic sur- 
veys. This methods requires a galaxy spectroscopic redshift sur- 
vey to reconstruct the large scale peculiar velocity and then weigh 
the galaxy density with this velocity to produce a weighted galaxy 
momentum map 0. By construction, this map is tightly correlated 
with the true kSZ signal and thus allows us to extract the kSZ sig- 
nal from the CMB maps through the cross cor relation measure - 
ment. Our work confirms the major finding of IHo et ah I ( l2009h . 
namely cross correlating the properly weighted galaxy momentum 
field with CMB fluctuations can significantly improve the kinetic 
SZ measurement. We estimate that BigBOSS-N plus PLANCK is 
able to measure the kSZ effect at 50cr confiden ce level, qualita- 
tively consistent with results of IHo et al] ( |2009|) . Given this S/N, 
our method is able to recover the redshift distribution of the kSZ 
effect over ~ 10 redshift bins. 

Thanks to the hydrodynamical simulations at hand, we are 
able to better model the true kSZ signal and thus reduce one source 
of error in the analysis. We are also able to quantify the impact of 
redshift distortion, by performing the reconstruction in both real 
and redshift space. We find that redshift distortion is a major fac- 
tor affecting the kSZ tomography. It degrades the performance, 
especially at small scales. However, we show that, its impact at 
k± 1/^/Mpc is moderate and thus the kSZ tomography can still 
work reasonably well to scales { =^ k±x ~ 3000. 

Comparing between the simulations with star formation, gas 
cooling and supernovae feedback turned on or off, we are able to 
quantify the impact of gastrophysics.We showed that the correla- 
tion coefficient r between the reconstructed and the true kSZ 
signal is insensitive to gastrophysics of star formation, feedback 
and gas cooling included in our hydrodynamic simulations. This 
behavior is of crucial importance for the theoretical interpretation 
of the kSZ tomography r esult. It enables a n approach similar to the 
thermal SZ tomography l lShao et alj|2009l) to circumvent the prob- 
lem of the potentially (and likely inevitably) large bias in © with 
respect to 0. This is an issue for further investigation. 

The proposed kSZ tomography has unique advantages over 
detecting the kinetic SZ eff^ect through auto-correlations. (1) By 
construction, the cross correlation measurement only picks up the 
kSZ component in CMB. Due to the lack of characteristic direc- 
tional dependence, primary CMB, the thermal SZ effect, dusty star 
forming galaxies and any noise sources of scalar nature do not bias 
the cross correlation measurement. It is thus a rather clean way 
to eliminate the otherwise overwhelming systematical errors in the 
kSZ measurement. (2) It recovers the redshift information of the 
kSZ effect. The redshift information is not only useful to better 
understand the evolution of missing baryons, but also useful to sep- 
arate the kSZ eff^ect after reionization from the one due to patchy 
reionization. It thus allows for better understanding of the reioniza- 
tion process. (3) The relatively higher S/N of galaxy surveys help 
to beat down the statistical errors in the kSZ measurement, as we 
have shown for the combination of PLANCK plus BigBOSS-N. 

We have carried out a concept study and shown the applicabil- 
ity of the kSZ tomography to real surveys. However, the analysis 
is simplified, with many issues left for future study. An incomplete 
list is as follows. 

(i) The galaxy stochasticity. It definitely degrades the kSZ tomog- 
raphy performance. However, studies shown that the stochastic- 
ity in the galaxy bias is likely at the level of 10% or smaller 
( iBonoh & Peij2009l ; lBaidauf et alj20l"ol) . So it is unlikely to com- 
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S iletely invalidate the kSZ tomography, consistent with IHo et alj 
2009f) in which the stochasticity exists in their galaxy mock cata- 
logue. Nonetheless, it is important to quantify the impact of galaxy 
stochastic bias. 

(ii) The smoothing scheme. The actual galaxy density distribution is 
irregular, due to the intrinsic nonlinearities, galaxy shot noise and 
survey irregularities. So an important step is to smooth the density 
field and stabilize the velocity reconstruction. We have done sim- 
ple study on the smoothing scheme in 33.3. II and found that this 
is indeed important. For example, we find that a spherical Gaus- 
sian s moothing does n ot work in redshift space as well as in real 
space. |Hoet2] (Eooi) investigated the Weiner filter to deal with 
galaxy shot noise. A comprehensive study is required to develop an 
optimal smoothing scheme. 

(iii) The cosmic variance. Our simulated box size is not large enough 
to quantify the cosmic variance, as the correlation length of velocity 
is ~ 100 Mpc . Simulations with larger box size are required to 
study this issue, along with the galaxy stochasticity and smoothing 
scheme. 

(iv) The robustness of modeling r, the cross coefficient r between the 
reconstructed and the true kSZ signal 0. It plays a central role 
in quantifying the performance of the kSZ tomography and in con- 
verting the measured cross correlation signal to the kSZ auto cor- 
relation power spectrum. We have measured the impact of gastro- 
physics on r and found that, gastrophysics of star formation, su- 
pemovae feedback and gas cooling would only change r by several 
percent on scales of interest. We shall check this result against more 
comprehensive investigations on gastrophysics. We have measured 
the dependence of r on the smoothing length and found a strong 
dependence. This implies that r in real surveys which have compli- 
cated masks and noise distribution could differ from what we have 
presented here. 

(v) Improvements on the forecast. The forecast can be improved by 
including other sources of contamination such as dusty star form- 
ing galaxies. It can also be extended to other combinations such 
as CMB experiments like ACT and SPT and many other galaxy 
spectroscopic redshift surveys. Dusty star forming galaxies do not 
have significant effect on the kSZ tomography based on PLANCK 
plus BigBOSS, since the applicable scale is ^ < 3000. At smaller 
scales approachable to ACT and SPT, dusty star forming galaxies 
overwhelm both the primary CMB and the thermal and kinetic SZ 
effect. So for these surveys, this source of error must be taken into 
account appropriately. 

(vi) Probing missing baryons with the kSZ tomography. Although it is 
one of major motivations and the most important applications of the 
kSZ tomography, we have not carried out any quantitative analysis 
on this aspect. Our hydro simulations have rich information on the 
IGM. However, the simulation box is not sufficiently large to ro- 
bustly quantify the contribution of the missing baryons to relevant 
statistics in the kSZ tomography. Furthermore, more comprehen- 
sive and robust treatments on gasdynamics and galaxy formation 
are demanded to reveal the connections between the kSZ effect and 
density and velocity of galaxies. Again, this requires more hydro- 
dynamical simulations, which will be investigated elsewhere. 



6 ACKNOWLEDGEMENT 

We thank S. Ho for his insightful questions and suggestions to im- 
prove the paper. This work is supported in part by the National 
Science Foundation of China (grant No. 10533030, 10673022, 
10821302, 10873027&10878001), the Knowledge Innovation Pro- 



gram of CAS (grant No. KJCX2-YW-T05 & KJCX3-SYW-N2), 
the 973 program grant No. 2007CB815401 & 2007CB8I5402 and 
the CAS/SAFEA International Partnership Program for Creative 
Research Teams. The simulations were done at Shanghai Super- 
computer Center by the supports of Chinese National 863 project 
(grant No.06AA0IA125). 



REFERENCES 

Atrio-Barandela, F, Mucket, J. P., Genova-Santos, R. 2008, ApJ, 
674, L61 

Baldauf, T, Smith, R. E., Seljak, U., Mandelbaum, R. 2010, PRD, 
81,063531 

Bemardeau, F, Colombi, S., Gaztafiaga, E., Scoccimarro, R. 
2002, Physics Reports, 367, 1 

Bonoh, S., Pen, U. L. 2009, MNRAS, 396, I6I0 

Bregman, J. N. 2007, Annual Review of Astronomy & Astro- 
physics, 45, 221 

Carlson, J., White, M., Padmanabhan, N. 2009, PRD, 80, 043531 

Cen, R., Ostriker, J. P 1999, ApJ, 514, 1 

Cen, R., Ostriker, J. P 2006, ApJ, 650, 560 

Chang, T.-C, Pen, U.-L., Peterson, J. B., McDonald, P 2008, 

Physical Review Letters, 100, 091303 
Dave, R., et al. 2001, ApJ, 552, 473 

DeDeo, S., Spergel, D. N., Trac, H. 2005, arXiv:astro-ph/0 51I060| 
Dore, O., Hennawi, J. F, Spergel, D. N. 2004, ApJ, 606, 46 
Fukugita, M., Peebles, P J. E. 2004, ApJ, 616, 643 
Fukugita, M., Peebles, P J. E. 2006, ApJ, 639, 590 
Genova-Santos, R., Atrio-Barandela, F, Mucket, J. P., Klar, J. S. 

2009, ApJ, 700, 447 
Hall, N. R., et al. 2010, ApJ, 718, 632 

Hemandez-Monteagudo, C, Trac, H., Verde, L., Jimenez, R. 

2006, ApJ, 652, LI 

Hernandez-Monteagudo, C, Ho, S. 2009, MNRAS, 398, 790 

Ho, S., Dedeo, S., Spergel, D. 2009, arXiv : 090372845 1 

Iliev, I. T, Pen, U.-L., Bond, J. R., Mellema, G., Shapiro, P R. 

2007, ApJ, 660, 933 

Jing, Y. P 2005, ApJ, 620, 559 

Jing, Y. P, Zhang, P, Lin, W. P, Gao, L., Springel, V. 2006, ApJ, 

640, LI 19 
Kaiser, N. 1987, MNRAS, 227, I 

Lin, W.P, Jing,Y.P, Mao,S., Gao,L., McCarthy,I.G., 2006, ApJ, 
651,636 

Ma, C.-P, Fry, J. N. 2002, Physical Review Letters, 88, 21 1301 
McQuinn, M., Furlanetto, S. R., Hemquist, L., Zahn, O., Zaldar- 

riaga, M. 2005, ApJ, 630, 643 
Nozawa, S., Itoh, N., Suda, Y, & Ohhata, Y. 2006, Nuovo Ci- 

mentoB Serie, 121,487 
Ostriker, J. P, Vishniac, E. T. 1986, ApJ, 306, L5I 
Peebles, P. J. E. 1980, The Large Scale Structure of the Universe, 

Princeton University Press, 1980. 
Pueblas, S., Scoccimarro, R. 2009, PRD, 80, 043504 
Schlegel, D. J., et al. 2009, arXiv:0904.0468] 
Scoccimarro, R. 2004, PRD, 70, 083007 
Shao, J., Zhang, P, Lin, W., Jing, Y. 2009,l arXiv:0903.5317l 
Springel, v., 2005, MNRAS, 364,1105 

Sunyaev, R. A., Zeldovich, Y. B. 1972, Comments on Astro- 
physics and Space Physics, 4, 173 
Sunyaev, R. A., Zeldovich, I. B. 1980, MNRAS, 190, 413 
Vishniac, E. T. 1987, ApJ, 322, 597 
White, M., Hernquist, L., Springel, V. 2002, ApJ, 579, 16 



14 Shao et al. 



2.5 
2.0 
"-a 1 .5 

1 .0 
0.5 



real space 


z = 


redshift space 


z = 0.53 






---z = 1.02 






_._..z = 2.08" 





0.1 1.0 
k/h/Mpc 



10.0 0.1 1.0 
k/h/Mpc 



10.0 



Figure Al. The reconstruction bias K in tlie real space (left panel) and K 
in the redshift space (right panel). In the real space, the reconstructed signal 
agree well with at high redshifts on large scales up to ki_ IftMpc"', 
while overestimates by almost 60% at recent epochs. However, in the 
redshift space, the resulting 0' powers ai'e boosted on the largest scales 
due to the Kaiser effect, while on smaller scales, FOG effect dominates and 
brings down the powers, the severest at z=0. 
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APPENDIX A: THE RECONSTRUCTION BIAS fig 

As we aim to reconstnact the kSZ map from the galaxy spectro- 
scopic survey, it's important to assess the reconstruction process. 
Quantitatively, we use r (Eq. |5) and (Eq. |6]l to figure out how 
well the reconstruction is. As explained in 32.21 of the main con- 
text, the S/N of the cross power spectrum (Eq. IT}, which is what 
we are most concerned about the future observations, is solely de- 
pendent on the determination of the cross correlation coefficient r, 
while a biased b^, i.e. the relative amplitude, does not influence 
the estimation of the signal to noise of the cross power spectrum. 
However, b^ is a direct measure of the reconstruction bias, and it 
does reflect how the reconstructed field suffers from nonlinearity 
and the redshift space distortion. So we here address the behavior 
of fog, both in the real space and the redshift space. 

The results of b'^ and b"^ are shown respectively in the left 
and right panel of Fig. lAll We caution the readers to pay attention 
only to the relative differences between real and redshift space, or 
between different redshifts. The absolute value of b^ could be mis- 
leading, due to the simplistic assumption of bg = 1 . 

In real space, tends to more over-estimate at lower red- 
shifts (&0 increases with redshift). In redshift space, it is the op- 
posite case. This is likely caused by the competition between the 
linear redshift distortion (the Kaiser effect) and the Finger of God 
effect. The Kaiser effect (Kaiser 198"^ enhances the clustering of 
matter along the line of sight and thus the derived velocity and 0'. 
This is the reason the redshift space ^0 at z = 2.08 is larger than the 
real space b^ a t the same epoch. O n the other hand, the Finger of 
God effect (see lScoccimarroir2004l for reviews of the redshift space 
distortion and the FOG effect) suppresses the redshift space matter 
overdensity and hence 0'. The Finger of God effect is amplified 
by the nonlinearities with respect to the Kaiser effect, resulting in 
decreasing b^ with redshift. It is also responsible for the smaller b^ 
in redshift space at z = comparing to the one in real space. 
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Figure Bl. The ratio of the power spectra of iJ-mode velocity and B-mode 
velocity. The results are the ratios of total B-mode powers and ZJ-mode pow- 
ers, different from average method of other figures, is compai'able to V£ 
beyond k ^ 2h Mpc^'at ; = 0, and this chai'acteristic scale moves to smaller 
scales at higher redshifts as velocity is subject to less nonlinearities at early 
epochs. 



APPENDIX B; 
VELOCITY 



THE INTERFERENCE OF THE "B"-MODE 



Our velocity estimator (Eq. [3), by construction, can only recover 
the irrotational part of the velocity (the so called "E"-mode veloc- 
ity. Since "B"-mode velocity (rotational part of the velocity, 
Vfi) also contributes to the kSZ effect, missing Vg degrades the re- 
construction and the kinetic SZ tomography. This is not a severe is- 
sue at large scales where Vg is negligible. However, at small scales 
multi-streaming develops and Vg begins to grow. Turbulence in the 
gas fluid, which can be amplified by gastrophysical processes like 
supernovae feedback, also contributes to Vg. To better understand 
the impact of Vg, we directly measure it from our simulations. 

The E-B modes can be obtained straightforwardly in Fourier 
space by 



V£(k) : 
Vg(k) : 



v,(k) • k , 

V„(k) - V£(k) . 



(Bl) 



The ratio of the two power spectra in our non-adiabatic run is shown 
in Fig. IB II As expected, the power spectrum of the "B"-mode ve- 
locity is much smaller than the "E"-mode velocity power spectrum 
on large scales. This is especially true at z > 1 and k < 1/i/Mpc. 
This is one of the major reasons for the reasonably good perfor- 
mance of the reconstruction at these scales and redshifts. The B- 
mode increases towards lower redshifts where the nonlinearities 
are stronger. At z = and k > l/iMpc"', it is already compara- 
ble to the E-mode. This behavior is largely responsible to poorer 
reconstruction alk > l/iMpc"', especially in redshift space (refer 
to 0X21 . 

Velocity measurement can be tricky in simulations, due to 
numerical artifacts, sampling bias in the velocity assignment and 
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Figure B2. The comparison of the gas velocity power spectrum ratio 
Pys IP\E non-adiabatic simulation (thick lines) and the adiabatic one 

(thin lines). Also shown are the ratios of dark matter velocity power spec- 
trum ratio (lines for adiabatic, open circles for non-adiabatic). The two sim- 
ulations have different softening length, with 9.75/i"' kpc (non-adiabatic) 
vs. 4.9/r ' kpc (adiabatic). The difference in softening length has virtually 
no impact on P^^ IP\e Amk matter So the significant changes in the 
gas velocity power spectrum between two simulations is very likely caused 
by differences in the gastrophysics other than difference in the softening 
length. We need to model "B"-mode velocity carefully in the existence of 
these gastrophysical processes, an important ingredient in the kSZ effect 
modeling. 



ferences, if any, between the two simulations are likely caused by 
difference in the softening length. We find the dark matter velocity 
has virtually no change in the two simulations. So the influence of 
softening length to the dark matter velocity is negligible. Although 
we are not able to directly quantify the influence of softening length 
to the gas velocity, this result suggests that the influence of soften- 
ing length should be also under control for the gas velocity. 

To the opposite, the gas velocity differs significantly between 
the two simulations, reflecting the significant influence of gastro- 
physics. At z = 2.08, the B-mode velocity of the non-adiabatic 
run is a factor of a few larger than that of the adiabatic run, likely 
caused by strong star formation and supernovae feedback at that 
epoch. Interestingly, this amplification in the B-mode velocity by 
gastrophysics at z = 2.08 is significantly weakened or even altered 
at z = 0. This is likely correlated with the fact that star formation 
rate decreases significantly from z = 2 to z = 0. 

The vorticity in the gas velocity is larger than that in the dark 
matter velocity at all redshifts, all scales and in both simulations. 
This difference is likely caused by gas viscosity and dissipation. 
It is an important issue to understand in order to improve the kSZ 
modeling. 



missing gastrophysics (e.g. IPueblas & Scoccimarroll2009t) . We do 
not attempt to perform a detailed study on these issues. Rather, we 
present a simple comparison between the two simulations avail- 
able. The two simulations share the same initial conditions and 
the same cosmology parameters. However, the gravitational soft- 
ening lengths are different. The one of the non-adiabatic run is 
9.75/;"' kpc, twice as large as the adiabatic one, 4.9/i"' kpc. With 
the two simulations, we can perform interesting observations on 
the following issues: 

• Difference between the velocity field of dark matter and that 
of gas. Understanding this difference helps improve models on the 
kSZ efl'ect. 

• The impact of gastrophysics on the gas velocity. Radiative and 
dynamical feedback can both cause turbulence in the gas fluid and 
thus affect the B-mode velocity. Comparing the gas velocity in the 
two simulations, we are able to get a handle on this issue. 

• The impact of the simulation softening length. The two simu- 
lations have different softening length. Comparing the dark matter 
velocity in the two simulations, we can estimate its impact. The 
reason that we do not use gas velocity to do the comparison is that 
the gas velocity is also affected by different gastrophysics in the 
two simulations. 

The ratios of two power spectra, Py^/Py^, are shown in Fig. 
IB2I Since dark matter velocity is insensitive to gastrophysics, dif- 



